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(57) ABSTRACT 

In the absence of a gyroscope, the attitude and attitude rate 
of a receiver can be determined using signals received by 
antennae on the receiver. Based on the signals received by 
the antennae, the Doppler difference between the signals is 
calculated. The Doppler difference may then be used to 
determine the attitude rate. With signals received from two 
signal sources by three antennae pairs, the three-dimensional 
attitude rate is determined. 

6 Claims, 2 Drawing Sheets 
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USING THE GLOBAL POSITIONING 
SATELLITE SYSTEM TO DETERMINE 
ATTITUDE RATES USING DOPPLER 
EFFECTS 

This application claims priority to U.S. Provisional 
Patent Application No. 60/225,030, filed Aug. 10, 2000. 

ORIGIN OF THE INVENTION 

The invention described herein was made by an employee 
of the United States Government, and may be manufactured 
and used by or for the Government for governmental pur- 
poses without the payment of any royalties. 

The present invention is drawn to determining the attitude 
and attitude rate of a receiver using signals from the Global 
Positioning System (GPS) and calculating the Doppler dif- 
ference between signals using multiple antennae on a 
receiver. 

BACKGROUND 

Global Positioning System devices typically use a 
phased-locked-loop to track the frequency of GPS satellites. 
The phase -locked-loop contains a “P” (proportional to error) 
and a “D” (proportional to error rate) term, the latter term 
being used to infer attitude rates. GPS attitude rate deter- 
mination using the “D” term of the phase-locked-loop is 
noisy. To get any useful information from it, the data has to 
go through a low-pass filter. This latter operation limits the 
high-frequency (high rates) that the system can handle. 

Alternatively, different devices, such as gyroscopes may 
be used. A variety of position-determining devices can be 
used in conjunction with numerical differentiation to yield 
satellite attitude rates such as star trackers and digital sun 
sensors. Gyroscopes are very accurate, but are yet another 
device for the satellite to carry. Also, numerically differen- 
tiating attitude -determination is inherently a very noisy 
operation. 

The following patents teach attitude or attitude rate deter- 
minations using various methods: U.S. Pat. Nos. 5,101,356; 
5,185,610; 5,525,998; 5,268,695; 5,543,804; 5,561,432; 

5,757,316; 5,844,521; 5,917,445; 5,933,110; 5,928,309; 

5,990,827; 6,101,430; 6,005,514; 6,078,284; 6,133,874; 

6,266,582; and “Development and Flight Demonstration of 
a GPS Receiver for Space” Dissertation by Edgar Glenn 
Lightsey. All of the above references are hereby incorpo- 
rated by reference. 

The present invention handles both high and low fre- 
quency attitude rate determination with high accuracy, for 
example, a zero crossing vector could be couple with a 
quartz-based counter, thereby yielding long period estima- 
tions to within microseconds). Using Doppler difference to 
determine attitude rate can replace the need for a gyroscope, 
thereby allowing cheaper attitude rate determination. 
Cheaper systems could use GPS-only for positioning, atti- 
tude and attitude rate. 

SUMMARY OF THE INVENTION 

The invention is drawn to determining the attitude and 
attitude rate of a receiver based on signals received by 
antennae on the receiver. Based on the signal received by the 
antennae, the Doppler difference between the signals is 
calculated. With signals received from two signal sources by 
three antennae pairs, the three-dimensional attitude rate is 
determined. 

BRIEF DESCRIPTION OF THE DRAWINGS 

FIG. 1 shows a system with one GPS satellite and two 
antennae. 
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FIG. 2 shows a system with one GPS satellite and four 
antennae. 

DETAILED DESCRIPTION 

5 Consider two antennae 1 and 2 on a receiver 3 and the 
signals they receive from one GPS satellite 10. If the two 
antennae are fixed relative to the GPS satellite, and ignoring 
effects such as multipath, they receive the same signal 
slightly phase shifted due to their separation. If the two 
antennae are rotated about the receiver centroid 4 between 
them, there will be a shift in frequency between the two 
signals due to the Doppler effect. One may extract the 
difference frequency by passing the two signals through a 
multiplier: 

sin (DjfX sin ^ 2 t=-Vi cos[(co 1 +co 2 )d+ 1 ^ cos[(( 0 1 -co 2 )d 

The sum of frequencies term is extremely high since GPS is 
broadcasting at gigahertz frequencies. Consequently, a low 
20 pass filter designed to remove it would easily allow mega- 
hertz rotation rates to be determined. The low frequency 
term is related to the rotation of the two antennae about their 
common centroid. This approach also requires cross- 
correlation with the pseudo-random noise (PRN) code of the 
25 GPS space vehicle under consideration, suitably phase 
shifted. This signal is already available on GPS receivers to 
support their tracking and acquisition functions. 

Equations 1-4 refer to FIG. 1, and define some basic 
variables and vectors which relate the antennae to a GPS 
30 satellite. 



r=U\\BL 21 \\ 

(1) 


v^v-oxV r ± 

(2) 

35 

v 2 =v-axV r 2 

(3) 


v 2 -v 1 =2(iiX(r 2 -r 1 ) 

(4) 


where 

40 r is the vector r 2 -r 1 

r 1 is the vector from the centroid of the receiver to antenna 

1 

r 2 is the vector from the centroid of the receiver to antenna 

2 

45 BL 21 is the distance between antenna 1 and antenna 2 
v is the velocity of the receiver relative to the GPS 
satellite 

co is the angular velocity vector for the satellite 
50 Vj is the velocity of antenna 1 relative to the GPS satellite 
v 2 is the velocity of antenna 2 relative to the GPS satellite 
a is the centroid of the receiver to GPS satellite normal 
vector 

Subtracting the two signals eliminates the common velocity, 
55 leaving a velocity dependent on the baseline vector between 
the two antennae and the angular velocity. 



f=f t (l-va/c) 

(5) 

60 

/rt=/*(l-((v-(aar r ± ) ■ a)lc ) 

(6) 


fr2=ftO--((V+MX > 2 ) -«)/ C ) 

(7) 


fn-fr2=f t (<*X(ri-r 2 ))/c 

(8) 


65 v is the velocity from the GPS satellite to the receiver 
a is the unit vector pointing from the receiver to the GPS 
satellite 
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f r is frequency of the signal received by the antenna 
f t is the transmitted frequency of the signal from the GPS 
satellite 

f rl is the modified frequency at r 1 
62 is the modified frequency at r 2 
c is the speed of light 

Equation 5 is the standard Doppler equation. Equation 8 
relates the three unknown variables in the 00 vector to the 
scalar difference frequency. Multiple satellites can be used to 10 
form a simultaneous set of three equations as is well-known, 
however, due to r, only the component of 00 normal to r can 
be determined. To get the complete 00 , three non-collinear 
antennae are required. Then three simultaneous, linear, inde- 
pendent equations are formed in various ways, using two ^ 
satellites and two pairs of antennae. 

FIG. 2 shows a four-antenna array. Atennae 21-23 must 
not be coplanar with antenna 24. The one satellite 20 and 
four-antenna array permit angular rate determination in the 
plane normal to the direction of propagation of its signal. In 2 o 
order to yield full three-dimensional determination of atti- 
tude rate two satellites (not collinear nor parallel with 
respect to the four antennae array) need to be used. 

The invention has been described using GPS signals as 
the signal source, however any other appropriate signal 2 s 
source may be used. Also, minimum means squared estimate 
and other known methods may be used to further reduce 
noise effects on attitude rate estimation. 

I claim: 

1. A method of determining the attitude rate of a receiver 30 
comprising: 

sending a first signal from a first signal source; 
receiving the first signal with a first antenna on a receiver; 
receiving the first signal with a second antenna on the 
receiver; 

sending a second signal from a second signal source; 
receiving the second signal with the first antenna; 
receiving the second signal with the second antenna; 
receiving the first signal with a third antenna on the 40 
receiver; 

receiving the second signal with the third antenna; 
determining the Doppler difference between the first 
signal and the second signal; and 45 

determining a component of the attitude rate based on the 
Doppler difference between the first signal and second 
signal, wherein the first, second and third antennae are 
non-collinear. 
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2. The method of claim 1, wherein the first signal source 
is a first Global Positioning System satellite. 

3. The method of claim 1, further comprising: 
sending a second signal from a second signal source; 
receiving the second signal with the first antenna; 
receiving the second signal with the second antenna; 
receiving the first signal with a third antenna on the 

receiver; 

receiving the second signal with the third antenna; 
receiving the first signal with a fourth antenna on the 
receiver; 

receiving the second signal with the fourth antenna; 
determining the Doppler difference from the received 
signals; 

determining the attitude rate based on the Doppler differ- 
ence between all of the received signals. 

4. The method of claim 3, wherein the second signal 
source is a second Global Positioning System satellite. 

5. The method of claim 3, wherein the first, second, third 
and fourth antennae are non-collinear and the fourth anten- 
nae is non-coplanar with the first second and third antennae. 

6 . A method of determining the attitude rate of a receiver 
comprising: 

sending a first signal from a first signal source; 
sending a second signal from a second signal source; 
receiving the first signal with a first antenna on a receiver; 
receiving the second signal with the first antenna; 
receiving the first signal with a second antenna on the 
receiver; 

receiving the second signal with the second antenna; 
receiving the first signal with a third antenna on the 
receiver; 

receiving the second signal with the third antenna; 
receiving the first signal with a fourth antenna on the 
receiver; 

receiving the second signal with the fourth antenna; 
determining the Doppler difference between the received 
signals; 

determining the three-dimensional attitude rate based on 
the Doppler difference between all of the received 
signals. 





